Protein engineering and the use of molecular modeling and simulation: the case of heterodimeric Fc engineering.

Computational and structure guided methods can make significant contributions to the development of solutions for difficult protein engineering problems, including the optimization of next generation of engineered antibodies. In this paper, we describe a contemporary industrial antibody engineering program, based on hypothesis-driven in silico protein optimization method. The foundational concepts and methods of computational protein engineering are discussed, and an example of a computational modeling and structure-guided protein engineering workflow is provided for the design of best-in-class heterodimeric Fc with high purity and favorable biophysical properties. We present the engineering rationale as well as structural and functional characterization data on these engineered designs.

Improving biophysical properties of a bispecific antibody scaffold to aid developability: quality by molecular design.

While the concept of Quality-by-Design is addressed at the upstream and downstream process development stages, we questioned whether there are advantages to addressing the issues of biologics quality early in the design of the molecule based on fundamental biophysical characterization, and thereby reduce complexities in the product development stages. Although limited number of bispecific therapeutics are in clinic, these developments have been plagued with difficulty in producing materials of sufficient quality and quantity for both preclinical and clinical studies. The engineered heterodimeric Fc is an industry-wide favorite scaffold for the design of bispecific protein therapeutics because of its structural, and potentially pharmacokinetic, similarity to the natural antibody. Development of molecules based on this concept, however, is challenged by the presence of potential homodimer contamination and stability loss relative to the natural Fc. We engineered a heterodimeric Fc with high heterodimeric specificity that also retains natural Fc-like biophysical properties, and demonstrate here that use of engineered Fc domains that mirror the natural system translates into an efficient and robust upstream stable cell line selection process as a first step toward a more developable therapeutic.

LC-MS characterization and purity assessment of a prototype bispecific antibody.

Bispecific IgG asymmetric (heterodimeric) antibodies offer enhanced therapeutic efficacy, but present unique challenges for drug development. These challenges are related to the proper assembly of heavy and light chains. Impurities such as symmetric (homodimeric) antibodies can arise with improper assembly. A new method to assess heterodimer purity of such bispecific antibody products is needed because traditional separation-based purity assays are unable to separate or quantify homodimer impurities. This paper presents a liquid chromatography-mass spectrometry (LC-MS)-based method for evaluating heterodimeric purity of a prototype asymmetric antibody containing two different heavy chains and two identical light chains. The heterodimer and independently expressed homodimeric standards were characterized by two complementary LC-MS techniques: Intact protein mass measurement of deglycosylated antibody and peptide map analyses. Intact protein mass analysis was used to check molecular integrity and composition. LC-MS(E) peptide mapping of Lys-C digests was used to verify protein sequences and characterize post-translational modifications, including C-terminal truncation species. Guided by the characterization results, a heterodimer purity assay was demonstrated by intact protein mass analysis of pure deglycosylated heterodimer spiked with each deglycosylated homodimeric standard. The assay was capable of detecting low levels (2%) of spiked homodimers in conjunction with co-eluting half antibodies and multiple mass species present in the homodimer standards and providing relative purity differences between samples. Detection of minor homodimer and half-antibody C-terminal truncation species at levels as low as 0.6% demonstrates the sensitivity of the method. This method is suitable for purity assessment of heterodimer samples during process and purification development of bispecific antibodies, e.g., clone selection.

Studies of base pair sequence effects on DNA solvation based on all-atom molecular dynamics simulations.

Detailed analyses of the sequence-dependent solvation and ion atmosphere of DNA are presented based on molecular dynamics (MD) simulations on all the 136 unique tetranucleotide steps obtained by the ABC consortium using the AMBER suite of programs. Significant sequence effects on solvation and ion localization were observed in these simulations. The results were compared to essentially all known experimental data on the subject. Proximity analysis was employed to highlight the sequence dependent differences in solvation and ion localization properties in the grooves of DNA. Comparison of the MD-calculated DNA structure with canonical A- and B-forms supports the idea that the G/C-rich sequences are closer to canonical A- than B-form structures, while the reverse is true for the poly A sequences, with the exception of the alternating ATAT sequence. Analysis of hydration density maps reveals that the flexibility of solute molecule has a significant effect on the nature of observed hydration. Energetic analysis of solute-solvent interactions based on proximity analysis of solvent reveals that the GC or CG base pairs interact more strongly with water molecules in the minor groove of DNA that the AT or TA base pairs, while the interactions of the AT or TA pairs in the major groove are stronger than those of the GC or CG pairs. Computation of solvent-accessible surface area of the nucleotide units in the simulated trajectories reveals that the similarity with results derived from analysis of a database of crystallographic structures is excellent. The MD trajectories tend to follow Manning's counterion condensation theory, presenting a region of condensed counterions within a radius of about 17 A from the DNA surface independent of sequence. The GC and CG pairs tend to associate with cations in the major groove of the DNA structure to a greater extent than the AT and TA pairs. Cation association is more frequent in the minor groove of AT than the GC pairs. In general, the observed water and ion atmosphere around the DNA sequences is the MD simulation is in good agreement with experimental observations.

Human subtilase SKI-1/S1P is a master regulator of the HCV Lifecycle and a potential host cell target for developing indirect-acting antiviral agents.

HCV infection is a major risk factor for liver cancer and liver transplantation worldwide. Overstimulation of host lipid metabolism in the liver by HCV-encoded proteins during viral infection creates a favorable environment for virus propagation and pathogenesis. In this study, we hypothesize that targeting cellular enzymes acting as master regulators of lipid homeostasis could represent a powerful approach to developing a novel class of broad-spectrum antivirals against infection associated with human Flaviviridae viruses such as hepatitis C virus (HCV), whose assembly and pathogenesis depend on interaction with lipid droplets (LDs). One such master regulator of cholesterol metabolic pathways is the host subtilisin/kexin-isozyme-1 (SKI-1)--or site-1 protease (S1P). SKI-1/S1P plays a critical role in the proteolytic activation of sterol regulatory element binding proteins (SREBPs), which control expression of the key enzymes of cholesterol and fatty-acid biosynthesis. Here we report the development of a SKI-1/S1P-specific protein-based inhibitor and its application to blocking the SREBP signaling cascade. We demonstrate that SKI-1/S1P inhibition effectively blocks HCV from establishing infection in hepatoma cells. The inhibitory mechanism is associated with a dramatic reduction in the abundance of neutral lipids, LDs, and the LD marker: adipose differentiation-related protein (ADRP)/perilipin 2. Reduction of LD formation inhibits virus assembly from infected cells. Importantly, we confirm that SKI-1/S1P is a key host factor for HCV infection by using a specific active, site-directed, small-molecule inhibitor of SKI-1/S1P: PF-429242. Our studies identify SKI-1/S1P as both a novel regulator of the HCV lifecycle and as a potential host-directed therapeutic target against HCV infection and liver steatosis. With identification of an increasing number of human viruses that use host LDs for infection, our results suggest that SKI-1/S1P inhibitors may allow development of novel broad-spectrum biopharmaceuticals that could lead to novel indirect-acting antiviral options with the current standard of care.

Spectroscopic and molecular dynamics evidence for a sequential mechanism for the A-to-B transition in DNA.

The A-to-B form transition has been examined in three DNA duplexes, d(CGCGAATTCGCG)(2), d(CGCGAATTGCGC), and d(CGCAAATTTCGC), using circular dichroism spectroscopy, ultraviolet resonance Raman (UVRR) spectroscopy, and molecular dynamics (MD) simulation. Circular dichroism spectra confirm that these molecules adopt the A form under conditions of reduced water activity. UVRR results, obtained under similar conditions, suggest that the transition involves a series of intermediate forms between A and B. Cooperative and distinct transitions were observed for the bases and the sugars. Independent MD simulations on d(CGCGAATTCGCG)(2) show a spontaneous change from the A to B form in aqueous solution and describe a kinetic model that agrees well with UVRR results. Based on these observations, we predict that the mechanism of the transition involves a series of A/B hybrid forms and is sequential in nature, similar to previous crystallographic studies of derivatized duplexes. A simulation in which waters were restrained in the major groove of B DNA shows a rapid, spontaneous change from B to A at reduced water activity. These results indicate that a quasiergodic sampling of the solvent distribution may be a problem in going from B to A at reduced water activity in the course of an MD simulation.

Prokaryotic gene finding based on physicochemical characteristics of codons calculated from molecular dynamics simulations.

An ab initio model for gene prediction in prokaryotic genomes is proposed based on physicochemical characteristics of codons calculated from molecular dynamics (MD) simulations. The model requires a specification of three calculated quantities for each codon: the double-helical trinucleotide base pairing energy, the base pair stacking energy, and an index of the propensity of a codon for protein-nucleic acid interactions. The base pairing and stacking energies for each codon are obtained from recently reported MD simulations on all unique tetranucleotide steps, and the third parameter is assigned based on the conjugate rule previously proposed to account for the wobble hypothesis with respect to degeneracies in the genetic code. The third interaction propensity parameter values correlate well with ab initio MD calculated solvation energies and flexibility of codon sequences as well as codon usage in genes and amino acid composition frequencies in approximately 175,000 protein sequences in the Swissprot database. Assignment of these three parameters for each codon enables the calculation of the magnitude and orientation of a cumulative three-dimensional vector for a DNA sequence of any length in each of the six genomic reading frames. Analysis of 372 genomes comprising approximately 350,000 genes shows that the orientations of the gene and nongene vectors are well differentiated and make a clear distinction feasible between genic and nongenic sequences at a level equivalent to or better than currently available knowledge-based models trained on the basis of empirical data, presenting a strong support for the possibility of a unique and useful physicochemical characterization of DNA sequences from codons to genomes.

CTL recognition of a protective immunodominant influenza A virus nucleoprotein epitope utilizes a highly restricted Vbeta but diverse Valpha repertoire: functional and structural implications.

To investigate protective immunity conferred by CTL against viral pathogens, we have analyzed CD8(+) T cell responses to the immunodominant nucleoprotein epitope (NP(366-374)) of influenza A virus in B6 mice during primary and secondary infections in vivo. Unlike the highly biased TCR Vbeta repertoire, the associated Valpha repertoire specific for the NP(366-374)/D(b) ligand is quite diverse. Nonetheless, certain public and conserved CDR3alpha clonotypes with distinct molecular signatures were identified. Pairing of public Valpha and Vbeta domains creates an alphabeta TCR heterodimer that binds efficiently to the NP(366-374)/D(b) ligand and stimulates T cell activation. In contrast, private TCRs, each comprising a distinct alpha chain paired with the same public beta chain, interact very differently. Molecular dynamics simulation reveals that the conformation and mobility of the shared Vbeta CDR loops are governed largely by the associated Valpha domains. These results provide insight into molecular principles regarding public versus private TCRs linked to immune surveillance after infection with influenza A virus.

Root mean square deviation probability analysis of molecular dynamics trajectories on DNA.

The comparison and detection of the commonalities and differences in multiple structural ensembles is an important step in the use of molecular simulations to gain insight into the conformation and dynamics of complex biomacromolecules. While the average structure is often employed as the representative of an ensemble of structures in such comparisons, dynamic molecular systems with multiple conformational substates call for a more accurate representation that captures the complete dynamical range of the ensemble. We present a probability analysis procedure based on the root-mean-square differences among the structural ensembles that efficiently and accurately performs the relevant comparison.

Structural bioinformatics of DNA: a web-based tool for the analysis of molecular dynamics results and structure prediction.

UNLABELLED: We report here the release of a web-based tool (MDDNA) to study and model the fine structural details of DNA on the basis of data extracted from a set of molecular dynamics (MD) trajectories of DNA sequences involving all the unique tetranucleotides. The dynamic web interface can be employed to analyze the first neighbor sequence context effects on the 10 unique dinucleotide steps of DNA. Functionality is included to build all atom models of any user-defined sequence based on the MD results. The backend of this interface is a relational database storing the conformational details of DNA obtained in 39 different MD simulation trajectories comprising all the 136 unique tetranucleotide steps. Examples of the use of this data to predict DNA structures are included. AVAILABILITY: http://humphry.chem.wesleyan.edu:8080/MDDNA. SUPPLEMENTARY INFORMATION: Supplementary data including color figures are available at Bioinformatics online.

Can free energy calculations be fast and accurate at the same time? Binding of low-affinity, non-peptide inhibitors to the SH2 domain of the src protein.

The usefulness of free-energy calculations in non-academic environments, in general, and in the pharmaceutical industry, in particular, is a long-time debated issue, often considered from the angle of cost/performance criteria. In the context of the rational drug design of low-affinity, non-peptide inhibitors to the SH2 domain of the (pp60)src tyrosine kinase, the continuing difficulties encountered in an attempt to obtain accurate free-energy estimates are addressed. free-energy calculations can provide a convincing answer, assuming that two key-requirements are fulfilled: (i) thorough sampling of the configurational space is necessary to minimize the statistical error, hence raising the question: to which extent can we sacrifice the computational effort, yet without jeopardizing the precision of the free-energy calculation? (ii) the sensitivity of binding free-energies to the parameters utilized imposes an appropriate parametrization of the potential energy function, especially for non-peptide molecules that are usually poorly described by multipurpose macromolecular force fields. Employing the free-energy perturbation method, accurate ranking, within +/-0.7 kcal/mol, is obtained in the case of four non-peptide mimes of a sequence recognized by the (pp60)src SH2 domain.

Molecular dynamics simulations of the 136 unique tetranucleotide sequences of DNA oligonucleotides. II: sequence context effects on the dynamical structures of the 10 unique dinucleotide steps.

Molecular dynamics (MD) simulations including water and counterions on B-DNA oligomers containing all 136 unique tetranucleotide basepair steps are reported. The objective is to obtain the calculated dynamical structure for at least two copies of each case, use the results to examine issues with regard to convergence and dynamical stability of MD on DNA, and determine the significance of sequence context effects on all unique dinucleotide steps. This information is essential to understand sequence effects on DNA structure and has implications on diverse problems in the structural biology of DNA. Calculations were carried out on the 136 cases embedded in 39 DNA oligomers with repeating tetranucleotide sequences, capped on both ends by GC pairs and each having a total length of 15 nucleotide pairs. All simulations were carried out using a well-defined state-of-the-art MD protocol, the AMBER suite of programs, and the parm94 force field. In a previous article (Beveridge et al. 2004. Biophysical Journal. 87:3799-3813), the research design, details of the simulation protocol, and informatics issues were described. Preliminary results from 15 ns MD trajectories were presented for the d(CpG) step in all 10 unique sequence contexts. The results indicated the sequence context effects to be small for this step, but revealed that MD on DNA at this length of trajectory is subject to surprisingly persistent cooperative transitions of the sugar-phosphate backbone torsion angles alpha and gamma. In this article, we report detailed analysis of the entire trajectory database and occurrence of various conformational substates and its impact on studies of context effects. The analysis reveals a possible direct correspondence between the sequence-dependent dynamical tendencies of DNA structure and the tendency to undergo transitions that "trap" them in nonstandard conformational substates. The difference in mean of the observed basepair step helicoidal parameter distribution with different flanking sequence sometimes differs by as much as one standard deviation, indicating that the extent of sequence effects could be significant. The observations reveal that the impact of a flexible dinucleotide such as CpG could extend beyond the immediate basepair neighbors. The results in general provide new insight into MD on DNA and the sequence-dependent dynamical structural characteristics of DNA.

Induced fit and the entropy of structural adaptation in the complexation of CAP and lambda-repressor with cognate DNA sequences.

Molecular dynamics (MD) simulations of 5 ns on protein-DNA complexes of catabolite-activator protein (CAP), lambda-repressor, and their corresponding uncomplexed protein and DNA, are reported. These cases represent two extremes of DNA bending, with CAP DNA bent severely and the lambda-operator nearly straight when complexed with protein. The calculations were performed using the AMBER suite of programs and the parm94 force field, validated for these studies by good agreement with experimental nuclear magnetic resonance data on DNA. An explicit computational model of structural adaptation and computation of the quasiharmonic entropy of association were obtained from the MD. The results indicate that, with respect to canonical B-form DNA, the extreme bending of the DNA in the complex with CAP is approximately 60% protein-induced and 40% intrinsic to the sequence-dependent structure of the free oligomer. The DNA in the complex is an energetically strained form, and the MD results are consistent with a conformational-capture mechanism. The calculated quasiharmonic entropy change accounts for the entropy difference between the two cases. The calculated entropy was decomposed into contributions from protein adaptation, DNA adaptation, and protein-DNA structural correlations. The origin of the entropy difference between CAP and lambda-repressor complexation arises more from the additional protein adaptation in the case of lambda, than to DNA bending and entropy contribution from DNA bending. The entropy arising from protein DNA cross-correlations, a contribution not previously discussed, is surprisingly large.

Structure and axis curvature in two dA6 x dT6 DNA oligonucleotides: comparison of molecular dynamics simulations with results from crystallography and NMR spectroscopy.

Molecular dynamics (MD) simulations have been performed on the A6 containing DNA dodecamers d(GGCAAAAAACGG) solved by NMR and d(CGCAAAAAAGCG) solved by crystallography. The experimental structures differ in the direction of axis bending and in other small but important aspects relevant to the DNA curvature problem. Five nanosecond MD simulations of each sequence have been performed, beginning with both the NMR and crystal forms as well as canonical B-form DNA. The results show that all simulations converge to a common form in close proximity to the observed NMR structure, indicating that the structure obtained in the crystal is likely a strained form due to packing effects. A-tracts in the MD model are essentially straight. The origin of axis curvature is found at pyrimidine-purine steps in the flanking sequences.

Molecular dynamics simulations of the 136 unique tetranucleotide sequences of DNA oligonucleotides. I. Research design and results on d(CpG) steps.

We describe herein a computationally intensive project aimed at carrying out molecular dynamics (MD) simulations including water and counterions on B-DNA oligomers containing all 136 unique tetranucleotide base sequences. This initiative was undertaken by an international collaborative effort involving nine research groups, the "Ascona B-DNA Consortium" (ABC). Calculations were carried out on the 136 cases imbedded in 39 DNA oligomers with repeating tetranucleotide sequences, capped on both ends by GC pairs and each having a total length of 15 nucleotide pairs. All MD simulations were carried out using a well-defined protocol, the AMBER suite of programs, and the parm94 force field. Phase I of the ABC project involves a total of approximately 0.6 mus of simulation for systems containing approximately 24,000 atoms. The resulting trajectories involve 600,000 coordinate sets and represent approximately 400 gigabytes of data. In this article, the research design, details of the simulation protocol, informatics issues, and the organization of the results into a web-accessible database are described. Preliminary results from 15-ns MD trajectories are presented for the d(CpG) step in its 10 unique sequence contexts, and issues of stability and convergence, the extent of quasiergodic problems, and the possibility of long-lived conformational substates are discussed.

Molecular dynamics simulations of DNA curvature and flexibility: helix phasing and premelting.

Recent studies of DNA axis curvature and flexibility based on molecular dynamics (MD) simulations on DNA are reviewed. The MD simulations are on DNA sequences up to 25 base pairs in length, including explicit consideration of counterions and waters in the computational model. MD studies are described for ApA steps, A-tracts, for sequences of A-tracts with helix phasing. In MD modeling, ApA steps and A-tracts in aqueous solution are essentially straight, relatively rigid, and exhibit the characteristic features associated with the B'-form of DNA. The results of MD modeling of A-tract oligonucleotides are validated by close accord with corresponding crystal structure results and nuclear magnetic resonance (NMR) nuclear Overhauser effect (NOE) and residual dipolar coupling (RDC) structures of d(CGCGAATTCGCG) and d(GGCAAAAAACGG). MD simulation successfully accounts for enhanced axis curvature in a set of three sequences with phased A-tracts studied to date. The primary origin of the axis curvature in the MD model is found at those pyrimidine/purine YpR "flexible hinge points" in a high roll, open hinge conformational substate. In the MD model of axis curvature in a DNA sequence with both phased A-tracts and YpR steps, the A-tracts appear to act as positioning elements that make the helix phasing more precise, and key YpR steps in the open hinge state serve as curvature elements. Our simulations on a phased A-tract sequence as a function of temperature show that the MD simulations exhibit a premelting transition in close accord with experiment, and predict that the mechanism involves a B'-to-B transition within A-tracts coupled with the prediction of a transition in key YpR steps from the high roll, open hinge, to a low roll, closed hinge substate. Diverse experimental observations on DNA curvature phenomena are examined in light of the MD model with no serious discrepancies. The collected MD results provide independent support for the "non-A-tract model" of DNA curvature. The "junction model" is indicated to be a special case of the non-A-tract model when there is a Y base at the 5' end of an A-tract. In accord with crystallography, the "ApA wedge model" is not supported by MD.